Late time anisotropy as an imprint of cosmological backreaction 
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Backreaction effects of the large scale structure on the background dynamics have been claimed 
to lead to a renormalization of the background dynamics that may account for the late time acceler- 
ation of the cosmic expansion. This article emphasizes that generically the averaged flow is locally 
anisotropic, a property that can be related to observation. Focusing on perturbation theory, the 
spatially averaged shear, that characterizes the anisotropy of the flow, is computed. It is shown that 
this shear arising from backreaction differs from a homogeneous shear: its time evolution is different 
and its amplitude is completely determined by the cosmological parameters and the matter power 
spectrum. It ranges within (2-37)% at a redshift of order 0.5 so that the isotropy of the Hubble flow 
may allow to constrain the backreaction approach to dark energy. 
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I. INTRODUCTION 

In the standard cosmological framework [TJ, one as- 
sumes that on large scales, the universe is well described 
by a spatially homogeneous and isotropic spacetime, at 
least at the background level, so that its dynamics is 
obtained from the Einstein equations for a Friedmann- 
Lemaitre (FL) metric. Structure formation is then de- 
scribed using perturbation theory, which has proven to be 
successful to understand the existing observations from 
cosmic microwave background anisotropics to the growth 
rate of the large scale structure. 

While most of the observations are compatible with 
the assumption of a spatially homogeneous and isotropic 
universe on large scales, they exhibit a clumpy distribu- 
tion of matter on small scales and at late time. This 
is at the heart of a lively debate concerning the magni- 
tude of the backreaction of the large scale structure on 
the background dynamics. While its magnitude depends 
on the averaging procedure [5] (see also Ref. [I]) and 
on the actual small scale geometry of the universe, it 
has been mostly estimated using perturbation theory at 
linear and second order either in synchronous gauge [SJ- 
[7] or in Newtonian gauge [7HT2] (see e.g. Ref. [13] for 
reviews). One important conclusion is the existence of 
ultraviolet divergences (see, for example, Refs. [HI ITS]) 
when particular observables are calculated, with no con- 
vincing regularization schemes proposed so far. This lets 
open the question of the magnitude of the backreaction 
and its ability to explain the late time acceleration of the 
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cosmic expansion. This is a rather controversial subject 
and, according to some authors [51 |S1 Qj5] the backreac- 
tion of present inhomogeneities might explain, by itself, 
cosmic acceleration while, according to others [T7]> the 
effect of inhomogeneities is totally negligible. 

This article lies on a very simple, and almost trivial, 
remark: the spatially averaged flow has no reason to be 
isotropic. While obvious, this fact has been hidden in the 
generally adopted backreaction procedure [2 , in which 
the averaged dynamics is presented in a form mimicking 
the Friedmann dynamics so that attention has been fo- 
cused on the volume averaged expansion factor and on 
the backreaction terms (see below for definitions). Late 
time growth of a spatial anisotropy may thus be a spe- 
cific signature of backreaction. The goal of this work is 
thus to estimate the level of anisotropy of the averaged 
flow and address the following questions: (1) Is the shear, 
which is related to the anisotropy of the Hubble flow (and 
thus observable), divergent in the ultraviolet? (2) How 
is the shear magnitude related to the backreaction mag- 
nitude? (3) Can a bound on the shear on cosmological 
scales allow us to set a bound on the backreaction? 

To that purpose, we start by recalling the averaging 
procedure in § [II] and propose to split the kinematical 
backreaction term into a term describing the anisotropy 
of the flow and a genuine backreaction term. Under that 
form the averaged dynamics of any spatially homoge- 
neous (but not necessarily isotropic) flow on a surface of 
homogeneity remains unchanged and the genuine backre- 
action term strictly vanishes. Using perturbation theory, 
the shear is computed at lowest order in perturbation 
theory in § III and then computed explicitly in § [V] after 
having discussed the choice of the matter power spec- 
trum in § [IV] We have actually used two power spectra, 
one motivated by observations (with a simple analytical 
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form allowing for an exact integration) and another mo- 
tivated by theory, relying on the initial conditions from 
inflation and the transfer function of the standard ACDM 
model, but which requires numerical integration. As we 
shall show the scalar shear is divergent in the ultravio- 
let, which raises several questions that are summarized in 
§ |VI| Appendix A summarizes the main results of linear 
perturbation theory around a spatially Euclidean homo- 
geneous background spacetime. A detailed description 
of how to implement the shear average is given in Ap- 
pendix B, and the normalization of the power spectra is 
described in Appendix C. 



B. Buchert's formalism in a nutshell 

In practice, perturbations grow significantly only dur- 
ing the matter-dominated era, so that one can restrict 
the analysis to dust-filled universes, eventually with a 
cosmological constant. In such a situation, and as it was 
assumed in the original work by Buchert [2], one can 
pick up a foliation defined by a geodesic congruence n M 
which accidentally coincides with the four-velocity of 
comoving observers, i.e. = u^. 1 

The shear tensor can be expressed as 



fxu q h^u® 



II. AVERAGED COSMOLOGICAL DYNAMICS 



and it follows that the scalar shear takes the form 



A. Spacetime foliation 



(2.9) 



The driving idea, as first proposed in Ref . [2] , has been 
to rely on the 1+3 splitting of the universe [TS] associ- 
ated with a general reference timelike congruence n M that 
defines a foliation of spacetime. The choice of this con- 
gruence is often referred to as a choice of observers. The 
3-dimensional spacelike hypersurfaces normal to can 
then be defined by the equation S(x,t) — Sq — 0, with 
S(x,t) a scalar field and So a constant. Then 



(-d^Sd v S g^y/ 2 ' 



(2.1) 



which is normalized as n^n^ 



-1. This allows us to 



define h^ v , the projector on these hypersurfaces, as 

h^v = g^u + n^rty (2.2) 

which satisfies by construction h^ph^ = h^ v and h^n^ = 
0. One can then define the expansion 0, shear a^ v and 
vorticity u^ t „ of the flow as 



0^ = h«h^ a n p 



(2.3) 
(2.4) 



The spatial average of any scalar quantity A on a do- 
main T> is then usually defined as [2] 



(A(t,x)h 



1 



\h\A(t,x)d i x 



(2.10) 



where Vx> is the volume of the domain, defined by the 
requirement that (l)x> = 1, and h is the determinant of 
the induced metric h^ u . Such a spatial average is asso- 
ciated with the general reference timelike congruence 
of Eq. ( |2.1[ ) if and only if the average is performed in the 
gauge where S(x,t) is homogeneous (see Refs. [2TJI |2"T] 
and Appendix B). 



The average (2.10) does not commute with the time 



derivative and gives origin to the so-called Buchert-Ehlers 
commutation rule |22j 

d t (A(t,x)) v - (d t A(t,x)) v = (QA) V - (G) V (A) V 

(2.11) 

where (0)x> is the spatial average of the expansion 
defined in Eq. ( |2.5[ ). We can then define an averaged scale 
factor a-p, related to the volume Vx> of the integration 
domain (normalized by a reference v olume scale Vp ) , by 
az> = (Vd/V-Dc) 1 / 3 . Considering Eq. (2.11), we then have 
the useful relation 



They are explicitly given by 



= V M n", 



= h a h@ 



ax> 



(2.12) 



(2.5) where it is important to emphasize that this scale factor 
describes the averaged dynamics of the flow and cannot 
be identified with a metric component of a FL spacetime. 
The bottom line of the averaging procedure proposed 

(2.6) 

in Ref. [2] is to allow for the spatial averaging of scalar 
quantities and to derive an averaged dynamical flow by 



UJ 



flu 



EE h«h?V [a 



0} 



(2.7) 



Indeed, the assumption of Eq. (|2.1| implies that the 



vorticity strictly vanishes, uiuv = 0. 



1 In general n M , which defines a general reference flow, and u^, 
which defines the four-velocity of the observers comoving with 
the matter, may be different (see Ref. (19) for details). 
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averaging the generalized Friedmann and Raychaudhuri 
equations to get 



a-D 



^{p)v + - i -\{Qv + {n) v ) (2.13) 
A On 

3 <^+3-f ^ 



where one has defined the kinematical backreaction term 
as 

2 



o pS ^((e 2 >x,-(e>|,)-2(a 2 )p 



(2.15) 



and (H)t> refers to the spatial average of the curvature of 
the 3-dimensional metric induced on the hypersurfaces. 
Since the average of the continuity equation, p + Qp = 0, 
leads to 

d t (p)v + (e)v(p)v = 0, (2.16) 

one easily concludes that (p)n oc a^, 3 . 

This formalism has been extended to allow for arbi- 
trary coordinate systems [T2j [23] and any choice of gauge 
and slicing [TO1|2"U] . 



C. Rewriting of the averaged equations 



In the form (2.13 2.14), the dynamics of the averaged 



scale factor mimics the Friedmann equations for a spa- 
tially homogeneous and isotropic FL spacetime. The 
dynamics of the averaged flow depends strongly on the 
backreaction terms Qn and (7t}n- These terms are re- 
lated each other by the integrability condition 



1 



d t (Qva v ) 



1 







(2.17) 



'v 



but there remains a freedom in the determination of Qx> ■ 
To this purpose, the value of Qn has been computed 
using perturbation theories in many studies (see, for ex- 
ample, Refs. [5HT2 1 [T4" l fT5] ) . and many alternative models 
have been p roposed (se e, for example, Ref. |2~4"]). 

The form ( 2.13 2.14 1 and the splitting of the backreac- 



tion terms is highly suggestive since they match a Fried- 
mann form. It is however intuitively clear that the aver- 
age of any spatially homogeneous flow should let this flow 
unchanged. In particular, we would like the kinematical 
backreaction term to strictly vanish when averaging any 
spatially homogeneous flow. To that purpose, we rewrite 
the system (2.13-2.14) in the strictly equivalent form 



an 



8irG . , A 
a-u I i 3 



&v 4-7tG , . A 

— = 5-(p)x> + -w 

an 3 3 



Qv 
3 ' 



defining a new kinematical backreaction term as 



Qv 



(0)1) 



(2.19) 



(2.20) 



Under this form, the averaged flow is described by a set 
of equations that mimics the flow of comoving observers 
in a spatially homegenous universe, but not necessarily 
isotropic spacetime. It follows that the backreaction term 
Qv is generated only by inhomogeneities and is thus re- 
lated to deviations from local homogeneity. 

Indeed, it is clear that for a FL universe, (<J 2 )v = 
while when averaging any Bianchi universe, (cr 2 )v = cr 2 
and Qv = 0. In both cases the set of equations (2.18 



2.19 1 has the homogeneous form, while the standard kine- 



matical backreaction Q-p is different from zero for the 
later case despite no genuine backreaction. With such a 
rewriting it is possible to compare the backreaction terms 
induced by the perturbations on a background flow with 
either FL or Bianchi symmetry |25j but more important 
to compare the magnitude of the anisotropy of the aver- 
aged flow with the true backreaction term. 

The averaging of spatially anisotropic cosmological 
models has been considered in Ref. [26]. Our defini- 
tion of Q-p, and the split of the kinematical backreaction 
term in Q-p and (<t 2 )t>, differs from their proposal. We 
emphasize that their definition relies on the use of (a) 2 ?, 
where a itself is not well-defined, contrary to a 2 . For 
instance, considering a FL background, the lowest order 
for a 2 will be quadratic in first order perturbations (see 
Appendix A and next section) which are stochastic fields 
and one needs to define properly the square root of these 
stochastic fields before ensemble average. Indeed, start- 
ing with a non-isotropic background, one can define a as 
\J (cr 2 )(°) at lowest order and then perform a perturba- 
tive expansion. Given this, we prefer to adopt a more 
straightforward and pedestrian rewriting of the standard 
equation. 



III. QUANTIFYING THE DEVIATION FROM 
ISOTROPY 



As explained above, the goal of this study is to quantify 
the anisotropy of the averaged flow. For that purpose, 
we start from a FL universe with perturbations in order 
to compute (<J 2 )v- In such an approach, one needs to 
be reminded that the perturbations are stochastic fields, 
usually with Gaussian initial conditions. It follows that 
the spatially averaged quantities are also stochastic quan- 
tities. 

Taking X as a function of the perturbations and (X)x> 
as its average on a spatial domain, then, from a theoreti- 
cal point of view, we only have access to the distribution 
of (A')-p, that is to (A)-p, which is the ensemble average 
of (X)x>. In this section, we compute the lowest contribu- 
tion to the shear in perturbation theory. The perturba- 
tion equations at linear order are summarized in § |IIIB| 
and the shear is explicitly computed in § III C 
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A. Background dynamics 

For simplicity, we consider a spatially Euclidean 
Friedmann-Lemaitre universe so that the late time dy- 
namics is dictated only by the pressureless matter and 
the cosmological constant with density parameters 



f2 m0 — 



87rGp m0 
3#o ' 



fl A Q = 



A 
3flo 



(3.1) 



that satisfy f2 m o + f^Ao = 1- The Friedmann equation 
then takes the usual form 



E{z) 



H\z) 

m 



^mo(l + ^) 3 + ^AO, 



the solution of which is 



a(t) oc sinh 2/3 ( ^n A0 H t 



(3.2) 



(3.3) 



The normalization to the Hubble constant today, Hq, 
implies that 



1/2 
AO 



„3/2 



2 v J ~ (l-fiAo) 1 /? - 

so that the redshift is given by 



l + z = 



sinh 2 / 3 (IVn^H t) 



(3.4) 



(3.5) 



B. Perturbation equations 

For our concern, it will be sufficient to assume that our 
Universe can be described by a spatially homogeneous 
and isotropic FL spacetime with perturbations. 

We decide to fix the gauge freedom by explicitly choos- 
ing the gauge in which the scalar S(x, t) is homogeneous. 
In this way we can use the simple definition of Eq. ( 2.10 ) 



for the spatial average of a general scalar (see Appendix 
B). As already mentioned, considering a ACDM model 
and the scalar S(x,t) which defines the reference hyper- 
surface comoving with the matter, namely choosing the 
matter comoving gauge in which one imposes that Tq = 
(with i — 1, . . . , 3), we choose observers which are also 
geodesic. We fix the remaining gauge freedom by im- 
posing that E — and Ei — (see Appendix A for the 
definitions). The metric thus takes the form 



ds 2 = a 2 [-d?? 2 + 2{d t B + B t )dr)dx 
+ (1 + 2C5, 



2E lj \ Ax'Ax 1 } (3.6) 



where B and C are scalars, Bi is a transverse vector 

{diB 1 = 0) and Ejj is a traceless and transverse tensor 

(diE lJ = = E\). The gauge invariant Bardeen poten- 
tials are then given by 

V = -C-HB, <5> = UB + B' (3.7) 



and the density contrast and peculiar velocity by 

5 C = 5 1 V = v = -B. (3.8) 

The perturbation equations reduce to the four Ein- 
stein equations and the Euler equation (again, we restrict 
ourselves to a dust plus cosmological constant universe, 
where the sound velocity c s is zero) 



V 2 * 



* = $ 



1 



2Af2j 



a p8' 



1 



(3.9) 
(3.10) 
(3.11) 



^" + mV + [2H' + H 2 ]y = (3.12) 
V +HV = -$, (3.13) 



with Mpj 2 = 8ttG. 



Scalar shear 



Following the definition (2.9) the scalar shear is by 



construction a second order quantity, 
(a 2 )^ = ( (7 2 ) (D = o. 



(3.14) 



As a consequence, it can be specified to its lowest or- 
der using only first order perturbation theory. Consider- 
ing the particular gauge choice specified in Eq. (3.6 1, we 



obtain (see Appendix A for the expression in a general 
gauge) 



2aV)( 2 > = /;,.,/r-< \^ri'.'r 



-2B td E 13 



(3.15) 

= diX 



where Bi = diB + Bi and we use the notation X , 
for any field X. 

At first order in perturbations, vector and tensor 
modes are negligible compared to scalar modes so that 
we can safely neglect them. The average on a domain 
T> is then obtained by inserting the expression (3.15) in 
Eq. (p0|, 



1 

Vd 



d 3 a; 
2a? 



1 



{V 2 Bf 



(3.16) 



Its ensemble average is easily obtained by shifting to 
Fourier space (see Appendix B for details) to get 



i 



d 3 k 



3a 2 (ry) J (2tt 



fc 4 |£? fc (77)| 2 (3.17) 



where, for any field X, Xk denotes its Fourier modes 
and X^X/,' = \X)~\ 2 S^(k + k'), because of statistical 
homogeneity and isotropy. The former integral can be 
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evaluated by the use of Eqs. (3.8), (3.10) and (3.13) so 
that the Fourier components of B and are related by 



a(rj) 



a(»/)¥fc(»/)d»/. 



(3.18) 



As already concluded in many works [Fl ll2lfT5] and shown 
here in the Appendix B, the ensemble average of a spatial 
average of any second order quantity on a domain T> does 
not depend on the size and shape of the domain at lowest 
order in perturbation theory. 

The anisotropy of the flow is characterized by the di- 
mensional quantity 



" ((0>p) 2 ' 
In perturbation theory it is given by 



3 H 2 



2 ((9 (2) )) P 

3 H 



so that at lowest order, it reduces to 
1 <(a 2 )C 2 % 



E = 



H 2 



(3.19) 



(3.20) 



(3.21) 



IV. CHOICE OF THE POWER-SPECTRUM 

To go further and compute the quantities derived in the 
previous section, we need to specify the power spectrum. 

We shall consider two power spectra. The first one re- 
lies on the inflationary prediction and the transfer func- 
tion of the standard ACDM model [57] , Our second 
choice is obtained from observation, and we decide to 
use the APM power spectrum [28] . Even if it is not the 
most up-to-date observational power spectrum, it gives 
a good description and has a particular simple form that 
allows us to perform the computations analytically. The 
two choices are complementary and will allow us to get a 
good estimate of the order of magnitude of the averaged 
shear. 

The spectrum Ps(k,r]) of the matter density contrast 
with Fourier components 6k(rj) — Sp m k(r))/p m (rj) is de- 
fined from the two-point correlation function 



and we also introduce the notation 

fc 3 



Vs 



2tt 2 



(4.1) 



(4.2) 



and use similar definitions for any field. 

The power-spectra of 8 and *S? are indeed related since 
Eq. (3.9) implies that 



1 

2^ 



1 



2M 2 



2 4 2 



Ps 



(4.3) 



A. First choice: theory based 

A first and natural choice is to use the power spectrum 
that is deduced from inflation and the transfer function 
computed in a ACDM model. Following Ref. [T2] the 
power spectrum takes the form 



(4.4) 



where A^j is the primordial power of the curvature per- 
turbation and is given by 



CMB 



71,-1 



with fccMB = 0.002Mpc , n s = 0.96 and 
A = 2.41 x 10~ 9 , 



(4.5) 



(4.6) 



so that it is CMB normalized [29] . 

For a ACDM model, the growing mode of the matter 
density perturbation is given in terms of a hypergeomet- 
ric function as [1 



D + {t)=a{t) 2 Fx 



1 , - ; — ; — sinh 2 | - x/^l^nHnt 
' 3 ' 6 ' \ 2 



so that the growth factor takes the form 



9{z) 



<h 



2^1 



1- 



1 11 



3' 6 ' (l + z) s 



The transfer function is parameterized as 

Lo(q) 



L (q) + C (q)q 2 



with 



Co(q) 
q being defined as 

Q 



L (q) = ln(2e + 1.8g) 
731 



14.2 



1 + 62.55' 



,(4.7) 
(4.8) 

(4.9) 

(4.10) 
(4.11) 

(4.12) 



13.41fc eq rj m0 /i 2 Mpc _1 

B. Second choice: observation based 

We can also use a form of the power spectrum moti- 
vated by observations and, for simplicity, we assume that 
its form is similar to the one obtained from the APM sur- 
vey [3H] 



P S (k) = A 



fc 



(fc c ) 



3/2 



(4.13) 



G 



with k c — (l/20)/iMpc . The normalization constant A 
is related to a\ 



1 



d 3 xS(x)W Rs (x) 



(4.14) 



where Wr 8 is a top-hat (spherical) window function with 
radius R$ = 8h~ 1 Mpc, a\ is typically of order unity to- 
day (see, e.g., Ref. [3D]) and it is equivalent to the vari- 
ance a\ defined in Appendix C. 

In the following, we decide to take, as value for ct 8 , 
the one obtained using as a starting point the theory 
based power-spectrum described above. Namely, in order 
for the two spectra to match, we compute erf with the 
inflationary spectrum, which is CMB normalized, and a 
set of cosmological parameters. We then impose that 
the value of a\ obtained from the APM spectrum is the 
same, which determines A. Assuming h = 0.7 and a 
standard ACDM model, we get A = 4.21 x 10 6 ct| Mpc 4 . 
The details are given in Appendix C. 

The power spectrum for \f, given by Eqs. (4.3| and 



(4.13), is almost constant in the infrared (k <C k c ) and 
behaves as k~ 3 in the ultraviolet. It thus has the same 
behaviour as the spectrum obtained from assuming stan- 
dard inflation in the case of a spectral index n s = 1 
in the infrared but differs in the ultraviolet. The two 
power spectra for the matter density contrast are com- 
pared in Fig. [T] for z = 0. As one can note, they are 
in extremely good agreement in the range of values of 
k from which we get the main contribution to the shear 
(see Eqs. (5.6-5.12)) and where the non-linear effects can 
still be neglected. 




i(r 4 



k(Mpc _1 ) 

FIG. 1: Comparison of the matter APM power spectrum 
Ps (dashed line) and the ACDM power spectum (solid line) 
assuming (fi m o,fiAo) = (0.27,0.73) and a spectral index 
n s = 0.96. 



V. SHEAR CALCULATION 

We are now in a position to actually compute the shear 
of the averaged flow (<r 2 )x> using these two power spec- 



tra. As already mentioned, the ensemble average of the 
spatially averaged shear loses reference to the averaging 
domain 2?, so that we can drop the dependence on T>. 

For a ACDM model, Eq. ( |3~9| implies that * oc D + /a 
so that 



9M 
9o 



(5.1) 



where go = g(z = 0). Considering Eq. (3.18) to evaluate 
B k one gets 

x 1 , s f Z / H Q g(z') 
Bk = - 1 + z)—=r* h (0) / dz> ° yK > . 

aoH a J zt H(z')(l + z') g 

The integral has to be performed numerically and the 
solution can be written as 

B k = -(1 + z)-*— ®k(0)I(z, Q m0 , f! A0 ) (5.3) 
cl 11q 

with 

i( z ,n m0 ,n A0 ) = fdz* f° 9{z>) (5.4) 

Jzt H{z'){l + z') g 

where Zi is an initial redshift deep in the matter era, the 
choice of which does not impact the final result. 



Going on we find that, starting from Eq. (3.17), the 
ensemble average of the shear is given by 



1 



1 



67r 2 a 2 \a a H a 
dfcfc 6 |* fc (0)| 2 



(i + z yi 2 (z,n m0 ,n A0 ) 



(5.5) 



A. APM power spectrum 



Using the expression of the power spectrum ( 4.3 ■ 4.13[ ) , 
the previous computation can be performed analytically 
to obtain 



^ = ^§|^(l + ^ 2 (^0,ttA0) 



dfc- 



[i + (k/k c y 



13/2 ■ 



(5.6) 



This integral is clearly convergent in the infrared but 
diverges in the ultraviolet. 2 It can be regularized in 



2 We should note that, on th e cont rary, the value of the expansion 
rate (ax>/ a x>) given in Eq. | |2.13| is not divergent in the ultravi- 
olet. This is a consequence of the fact that the divergent term in 
the shear cancels with a similar divergence in the genui ne ki nc- 
matical backreaction (which at leading order is Qjj ~ 2{cr 2 ) T> ). 
Therefore, this divergent behavior can imprint on the level of 
anisotropy of the averaged flow while not affecting the total back- 
reaction on {axi/a-v) if calculated with the standard Buchert 
procedure. 
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two different ways: either by introducing an UV cut- 
off fcuv, or with the method introduced in Ref. [12] 
in which |* fc (0)| 2 is replaced by W{kR s ) 2 \f k (0)\ 2 with 
W(kR s ) — e~ k R =/ 2 . R s is a smoothing scale (and we 
will assume that R s = fcuy). 

The result of the integration then takes the form 



3H 2 



47T 2 



n 2 m0 f(z)w(k vv ) 



(5.7) 



Ak 9 

where ^— f f2 ln0 is a constant but depends on the partic- 
ular set of cosmological parameters considered, the red- 
shift evolution is characterized by 



/(*) 



(1 + ^ 
E(z) 



i 2 {z,n m0l n A0 ) 



(5.8) 



and IV (fcuv) is a function depending on the regulariza- 
tion scheme and the cut-off (but not on the cosmological 
parameters). With a simple cut-off 



^%v) = -l+(l + §|)(l 



fc 2 

fc 2 



-1/2 



(5.9) 



while the regularization based on a Gaussian window 
function gives 



W(k 



uvj — 



s/tt k 



^ e fe c 2 /4v I 1 



k 2 

K UV- 



1-Erf 



(5.10) 



where Erf is the error function. 



for a window function. In both cases, the integration 
cannot be performed analytically and we have to rely on 
a numerical integration. Note, however, that the result 
will depend only mildly on the spectral index as long as 
n s ~ 1. 

These two window functions are compared on Fig. [2] 
considering both the APM and the inflationary power 
spectrums. In both cases, it can be concluded that 
the dependence on the choice of the window function, 
Eqs. (5.9) or ( 5.10 1 for the APM power spectrum and 
Eqs. (5.131 or (5.14) for the inflationary power spec- 
trum, is mild. As long as fcuv remains small compared 
to IMpc -1 , there is almost no difference. See also Fig. [3] 
below, where we have depicted, for the APM case, the 
time evolution of the shear for the two window func- 
tions assuming fcuv = 0.1 Mpc -1 , fcuv = 0.5 Mpc -1 and 
fc uv = lMpc^ 1 . 



>S 0.1 - 




(Mpc -1 ) 



B. ACDM power spectrum 

Let us now repeat this calculation with the power spec- 
trum ( 4.4 ) . The expression of the shear is now given by 



2 A k: / 5o xJ 





3H 2 25 alH* \g 
where we now have 



f(z)W(k vv ) (5.11) 



n,-l 



This expression is clearly convergent in the infrared while 
it diverges in the ultraviolet so that we have to regularize 
it. 

Proceeding as in the previous paragraph, the use of an 
UV cut-off leads to 



W 



1 r kvv dfcfc 3 / fc 



2-/o ^ c 



"CMB 



T 2 (fc), (5.13) 



1 TdfcfcV * ) ( a ^vr r , ( , !( ,. 14 , 

2 Jo k c V «CMB 



FIG. 2: Dependence of the window function W on the UV 
cut-off scales fcuv, assuming the APM (thin lines) or an infla- 
tionary (thick lines) power spectrum. The solid line refers to a 
sharp cut-off and the dashed line to a Gaussian regularization. 



C. Properties of the shear 

The two previous sections show that the choices of the 
cut-off scale and of the power spectrum affect the overall 
value of the shear. However, since, as seen from Fig. [T] 
the two spectra are very similar in shape, we do not have 
an important difference between the two choices. It is 
also important to realize that the redshift dependence of 
the shear is dictated by the function f(z), which is the 
same in the two cases. 

Let us first discuss the amplitude. The typical am- 
plitude of £ is almost unaffected by the choice of power 
spectrum. Typically the value of S computed with the 
inflationary spectrum will be from 3% to 9% smaller than 
the value computed with the APM spectrum when fcuv 
changes from 0.1 Mpc -1 to IMpc -1 . For that reason, 
we will now present just one set of results. Fig. [3] de- 
picts the evolution of E for different values of fcuv and it 



can easily be checked that the ratio between the overall 
amplitudes of two curves coincides with the ratio of the 
corresponding two window functions W(k\jv). Typically, 



the overall amplitude is multiplied by a factor of order 
9 or 19, respectively, when using /cyv = 0.5 Mpc -1 or 
k\jv = IMpc -1 instead of fc uv = 0.1 Mpc -1 . 




Let us now discuss the time evolution of E, which 
follows from the form of the function f(z) defined in 
Eq. (5.8). The evolution of I(z, O m o, Qao) at large red- 
shift (typically z > 1) can be obtained by assuming that 
g(z) is almost constant and g(z) ~ goo. In this case we 
have 



j i o o \ 1 5oo r) -i/2 1 \A + a o(l + zf -1 
3 go y/l + a {l + z) 3 + 1 



with ao — fl m o/Q\Q. It follows that 

(a 2 ) v Aki (9oo\ 2 w(h , 1 



(5.15) 



for z > 1. This is indeed consistent with our numerical 
results, as depicted on Fig. [3] 

In 2 = 0, the function / can be expanded as f(z) ~ 
f(0) + f'(0)z so that the function starts increasing only if 
/'(0) > 0, in which case E(z) will peak at a redshift z max - 
The condition for the existence of such a peak takes the 
form 



1(0, ^mOi ^Ao) 



> 



(5.16) 



since from Eq. (5.4), J'(0, O m o, fl/vo) = 1- This condition 
cannot be easily translated into a condition on the cos- 
mological parameters, but for a flat ACDM, one gets that 
there exists a peak if SIao > 0.44. Interestingly, the posi- 
tion of the peak depends on the value of the cosmological 
constant. Fig. [3] depicts the value z max of the position of 
the peak as a function of Oao This computation can 
be compared to the explicit form of Ti(z) allowing Qao 
to vary between 0.2 and 0.9; see Fig. [3J Note that one 
needs to compute A for each value of the cosmological 
parameters since the relation between cr| and the the- 
ory based power spectrum depends on the growth rate 



and fc G q- The amplitude of the peak decreases with 17ao 
which can be traced back to (1) the fact that the growth 
rate is smaller for larger A and (2) the fact that k eq de- 
creases with f^Ao- We will have a longer radiation era 
with more modes entering the Hubble radius during this 
era. As a consequence, they will be significantly damped 
with respect to those which become sub-Hubble during 
the matter era (see, for example, Ref. [T5]). 




FIG. 4: Value z max of the position of the peak of E as a 
function of SIao- Such a peak exists only if the condition ( 5.16 1 
is satisfied, that is for Qao > 0.44. 

The behavior of has to be compared with the 

behavior of a homogeneous shear, e.g., in a Bianchi I 
universe, for which it scales as a 2 oc (1 + z) 6 , so that in 
the matter era, it decays as Eeianchi ~ z 3 |31j while the 
spatially average anisotropy induced by the large scale 



structure grows as S 



to eventually peak and then 



slowly decay as S ~ z. Its typical amplitude is set by the 
value of the cosmological parameters and the normaliza- 
tion of the matter power spectrum, contrary to a Bianchi 
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FIG. 5: Dependence of £(2) on the value of JIao for an UV 
cut-off fcuv = 0.5 Mpc -1 . The position of the peak corre- 
sponds to the value of z max depicted in Fig. [4] JIao ranges 
from 0.2 to 0.9 by steps of 0.1 (with the addition of the curve 
correspondent to Q ao = 0.73). The dashed curves do not sat- 
isfy the condition (5.161 and have no maximum. The thick 
curve is the one which corresponds to Q,ao = 0.73. 
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FIG. 6: Amplitude of the peak, £(z max ) as a function of £Iao- 
From bottom to top fcuv = 0.1, 0.5, 1 Mpc^ 1 and we have used 



the two regularization schemes (solution ( 5.9 1 (solid lines) and 
solution (5.101 (dashed lines)). It decreases with Oao- 



universe for which it is a pure initial condition. 



VI. DISCUSSION 

This article relies on the simple remark that generically 
a spatially averaged flow has no reason to be isotropic. 
We have thus focused our analysis of the shear in order to 
quantify the expected deviation from isotropy that arises 
from the backreaction of large scale structure. 

For that purpose, we have used the averaging proce- 
dure designed in Ref. [2] but rewritten in a way that 
makes explicit the anisotropy of the averaged flow. 
Under this new form the kinematical backreaction term 
arises only from inhomogeneities and, e.g., the spatial 
average of a Bianchi dynamics remains unchanged, so 
that it is physically more sound (in particular, one 
can think of comparing the backreaction that arises 



from similar inhomogeneities with respect to different 
background homogeneous flows). In order to evaluate 
the order of magnitude of the scalar anisotropy S, it is 
sufficient to work at first order in perturbation theory. 

Concerning the questions that were raised in the intro- 
duction, we can give the following answers. 

First, (ex 2 ) diverges in the ultraviolet and it must be 
regul arize d by t he int roduction of a U V cu t-off fc uy; see 
Eqs. Kty and (|5.10b or Eqs. (|5.13|) and (I5.14I). Such 



UV divergence has a large impact on the final result and 
the same problem of UV divergent terms is also found in 
other approaches (e.g. the one based on the all-sky aver- 
age (monopole) of the redshift-distance relation, obtained 
from the Kristan and Sachs approach [35], as clearly 
shown in Ref. _15j). 3 The question regarding this cut- 
off is central in the backreaction debate. Indeed, since 
we have used linear perturbation theory, it is probably 
not safe to extrapolate at scales larger than k ~ 1 Mpc -1 
without taking into account the non-linear effects which 
turn on beyond such scales. This divergence is particu- 
lary problematic for the shear, since this is related to an 
anisotropy of the Hubble flow and can thus, in principle, 
be observed. We see that its amplitude can shift from 
2% to 37% when fcrjv varies from 0.1 to 1 Mpc -1 for the 
standard ACDM model (see Fig. 3), and it varies with 
the value of the cosmological constant (see Fig. 6) . 

Since the amplitude of some backreaction effects can 
also depend on /ctjv (see Ref. [H]), a constraint of the 
scalar shear can be translated into a constraint on the 
backreaction, hence possibly resolving the debate of the 
ability for these effects to explain the late time accel- 
eration of the universe. The average Friedmann equa- 
tion ( 2.18 ) can be rewritten as 1 = Q m + ft a + Q<j + + 
CIq. It is well-known (see, for example, Refs. [5j[T2]) that 
£Iq = fl a + Oq, which usually characterizes the backre- 
action, is not UV divergent, while both and J7q are. 
This arises from the fact that the divergence does not 
appear in (0)x> while it does in both (a)u and (0 2 )x>- 
This shows trivially that il a ~ — £Iq for large fcrjv so * na * 
the amplitude of the anisotropy and kinematical backre- 
action are related in this regime. 

Let us stress that there is no total agreement about the 
physical meaning of the Buchert procedure and its con- 
nection to observable quantities. For example, quanti- 
ties connected to spatial averages are not, strictly speak- 
ing, observable because we only have observational access 
to our past light-cone (see, for example, Ref. 15J). On 
the other hand, if we take seriously these equations we 
have an important observational impact via the averaged 
shear. 



3 However, note that the use of a well-defined procedure to average 
over the light-cone I33| , applied to the perturbative expansion up 
to the second order of the exact redshift-distance relation, gives 
a well-defined result with no ultraviolet divergences 1341 . 
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It is then important to realize that the anisotropy aris- 
ing from the averaged flow, as a specific signature of 
backreaction, may allow us to discriminate with dark en- 
ergy models |35j (unless dark energy has an anisotropic 
stress, as e.g. in Ref. [36 ) or modification of general 
relativity. We argue that this is a key quantity to dis- 
criminate these models as a possible explanation of the 
late time acceleration of the cosmic expansion. We have 
shown that its amplitude is fixed by the cosmological pa- 
rameters and the matter power spectrum, and not as an 
initial condition as e.g. for Bianchi universes. Its time 
evolution is also different from the homogeneous shear of 
a Bianchi universe since it increases with time to even- 
tually reach a maximum. In an effective way, this can 
indeed be rephrased as a dark energy anisotropic stress, 
the time evolution of which is determined by the growth 
of the structures. 

We have shown that for a standard ACDM background 
dynamics, £ peaks between 2 and 37% at a redshift of 
order 0.5. Any bound on the anisotropy of the Hub- 
ble flow can probably set constraints on backreaction 
effects. In particular, it shows explicitly that a cut-off 
scales larger than or of the order of 1 Mpc -1 induces too 
strong an anisotropy. Indeed, it does not make sense to 
derive the value of a cut-off from observation, but it em- 
phasizes, again, the necessity to design a proper regular- 
ization scheme for this mechanism to be fully predictive. 
An anisotropy with that amplitude can potentially be 
observable, in particular with weak lensing experiment 
such as Euclid [37] via the B- modes [55] . 

From a technical point of view, we have used two mat- 
ter power spectra. One is more realistic while the second 
has a more simple form. Interestingly, we have shown 
that they give similar results as long as k\jy remains 
small. This is thus a excellent tool that allows us to make 
analytical estimates and, for instance, to derive scalings 
and order of magnitude without numerical integrations 
of the spectrum. 

Our analysis, while restricted to simple power spectra 
and linear perturbation theory, points toward a specific 
signature of backreaction. This may offer the possibility 
to constrain this class of explanations for the dark energy. 

Acknowledgements: We thank Chris Clarkson, 
Obinnah Umeh and Carlo Schimd for discussions and 
the Action Specifique GRAM for its financial support. 
We also thank Thomas Buchert for his comments and 
insights. 



Appendix A: Linear peturbation theory 

We summarize here the main result of linear pertur- 
bation theory around a spatially Euclidean FL space- 
time. We use the notation of Ref. [1 . the standard 
scalar-vector-tensor decomposition of the metric compo- 



nent takes the form 
S^g 0Q = -2a 2 A, 
5 (1) 5i0 = a 2 B i = a 2 (d i B + B i 



5^ 9ij = a 2 [S l3 2C + 2d t d ] E + 2d (l E j) 



2E„. 



with 4 scalar degrees of freedom (A, B, C and E), 2 
transverse vectors (Bi and Ei with d l Bi = 0, d l Ei = 0) 
with 4 degrees of freedom, and a traceless and transverse 
tensor (Eij with d % Eij = = E\) with 2 degrees of free- 
dom. One can then define 6 gauge invariant degrees of 
freedom usually defined as 

* = -C -H{B- E'), (Al) 

$ = A + H(B -E') + (B-E'Y , (A2) 

¥ ee E 1 ' - & , (A3) 

E ij (A4) 

where the prime denotes the derivative with respect to 
conformal time and % = a' / a. 

Considering a matter sector described by a perfect fluid 
with stress-energy tensor 



Tftu = (p + P^u,, + Pg^ v , 



(A5) 



where the density and pressure can be splitted as 
p(x, rj) = p(rj) + Sp{x, rj) and P(x, rj) = P(rj) + SP(x, r/), 
and the velocity of the comoving observers is decomposed 
as it M = + 5u^ with u u u^ = — 1. It follows that 



u. 



a{-l~ A,v l + B l ) (A6) 



and we decompose Vi into scalar and a vector component 
according to 



diV + Vi . 



(A7) 



Some of the gauge invariant variables associated to the 
matter sector are then given by 



S c = 6 

V = v 
Vi = & 



P_ 
P 
E' 

-Bi 



(v + B) 



(A8) 

(A9) 
(A10) 



The scalar shear, given in Eq. (2.9), is by construction 
a second order quantity so that 

(a 2 )^ = (a 2 )^ - 0. (All) 

Thus, at the lowest order and for a general metric, we 
have 



2n(2) =: 



1 



+ 



2a 2 S' 2 
1 
2^2 
1 

1 



SS^B 1 ' 3 



1 



(V 2 8Sf 



1 



(V 2 6S)(d l B z ) 



a 2 S> 
2a 2 ' 



SSah ^ Ji, ,h ■'■> 



(A12) 
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\5 i3 ;S7 2 E + d l E 3 + djEi + E i3 , and 



where h_^i ■ : \ * i >>,i..,, <// 

diX for any field X 
It is clear that first order perturbation theory is suffi- 
cient to obtain the general expression for the shear up to 
second order, since second order perturbations will con- 
tribute only to third or fourth order to a 2 . 

Let us stress that if we neglect vector and tensor per- 
turbations at first order (considering their contribution 
negligible with respect to the scalar one) then there will 
be no vector or tensor contributions in (a 2 Y 2 \ In fact, 
genuine second order vector and tensor contributions, 
that are always sourced by first order scalar perturba- 
tions, are also not present in (er 2 )( 2 ) as we have seen. We 
can note that this is also a consequence of a general prop- 
erty valid for any scalar field at any order in perturbation 
theory. Namely, at a given order n the vectors B\ and 
E\ n \ and the tensor can appear in a scalar quan- 



tity only as d i Bl" ) , d'E^', <V, 7' //,'/' or h\ {n \ but since 
the vectors and the tensors are, respectively, transverse, 
and traceless and transverse to all order, these terms are 
always identically zero. 



manifestly gauge invariant form, which allows the use of 
different gauges independently of the choice of the space- 
like hypersurface on which the average is performed. For 
instance, the average of a scalar quantity A(t, x) on a 
hypersurface S(t, x) = So can be obtained from the four- 
dimensional integral 



I(A,fl) = [ d i x^/-g(t,x)A(t,x) 
Jn(t,x) 

d 4 xy/-g(t, x) A{t, x)W a {t, x) . (Bl) 



Here g is the determinant of the 4-metric g^ v and the 
window function is given by 



Wn(t, x) = n"V M 0[S(i, x) - S ]W n (t, x), (B2) 



Appendix B: Shear average 

Following Refs. p~9l I20| . the spatial average prescrip- 
tion described in Eq. (2.10) can be generalized to take a 



where Wn(t,x) defines the 3-dimensional domain D in- 
side the 3-dimensional hypersurface S(t,x) — Sq. Then, 
the average is simply given by [HI [20] 



( D (x)){S ,r } 



f Ssg d 3 x^\i(t ,x)\ D{t Q ,x)W n [h{t Q ,x),x)} 
JjT d3 Xy /Wt^\Wn[h{t Q ,xj^] ' 



(B3) 



r 



where 7 is the determinant of the 3-metric gij and the bar 
indicates that the quantities are evaluated in a new bar 
coordinate system in which the scalar S is homogeneous 
(i.e. S(t,x) = S(°'(t)). This corresponds to a redefi- 
nition of the time coordinate, t — > t with t = h(t,x). 
Furthermore, ta is the time t at which S^(t) takes the 
constant values S^. The suffix £s on the integral recalls 
that we are averaging a scalar quantity A on a section of 
the three-dimensional hypersurface Ss D , where the given 
scalar field S(t, x) takes the constant values So- 

On the other hand, for the aim of the paper, the sim- 



,(2) 



(a 2 f\ 



invariant up to this second order, i.e. (a 
Therefore there is no difference between the particular 
bar gauge and any other gauge for such a quantity. As 
a final result, we obtain that, up to second order in 
perturbation theory, the average of the shear over a 3- 
dimensional domain T> can be simply written as 



1 



v Jv 



d 3 x(a 2 )^(t 0l x). 



(B4) 



Considering now the gauge introduced in Eq. (|3.6[) of 



III B using Eq. (3.15 1 and neglecting tensor and vector 



pier definition given by Eq. ( |2.10[ ) is sufficient to average contributions, we end up with 



the scalar shear without running into problems related to 
the gauge choice. In fact, a double simplification appears 
as a consequence of the fact that the scalar shear is non- 
vanishing only from the second order on (see Eqs.(3.14) 
and (3.15)). First, performing the computations only up 

to second order, we can take 



1 



2 u \2 



(B5) 



71 = Vl7 (0) l = a(t) 3 , 
and the particular form of the perturbations of the spa- 
tial metric is irrelevant. Second, the shear will be gauge 



Let us now consider the Fourier expansion of the first 
order perturbation B as 



B(x,r]) 



d 3 k 

(2tt) 3 /- 



e lkx B k ( V ). 



(B6) 
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The integral becomes 



lk' 2 k" 2 ) B k ,B k „ e « k '+ k "> x 



(B7) 



If, as predicted from inflation, the cosmological perturba- 
tions are related to the primordial perturbations which 
enjoy Gaussian statistics with zero mean pQ, the aver- 
age can be computed by taking an ensemble average over 
many domains 39J . We denote this additional average 
by an over-bar. As a consequence, if B is statistically 
homogeneous then 

B{k)B(k') = \B k \8 {3} (k + k'). (B8) 

The ensemble average of (<J 2 )v is thus given by 



W 2 h 



1 r ^ 1 1 

d^x- 



V v Jv ~ 2a 2 (2tt) 3 

1 1 
3a 2 (2tt) 3 



d 3 fc'' 



\k'\ A \B k ,\ 2 



d a k'\k'\*\B k 



(B9) 



and it will be independent from the domain of integration 
T>, whatever is the window function which defines such 
domain of integration. This is a general property for the 
ensemble average of the spatial average of a second order 
contribution. However, it is worth underlining that the 
result (B9| holds only because the higher contributions 



arising from the ensemble average, and involving also the 
induced metric (which is also present in the definition of 
Vp), will be at least fourth order in perturbation theory. 



Appendix C: Normalization of the matter power 
spectra 

1. Definition of <j\ 

For any density field 8, that can be decomposed in 
Fourier modes as 



8(x,t) 



d 3 k 

(2tt) 3 / 2 



8 k e l 



(CI) 



and any window function Wr(x) of typical width R, the 
smoothed density field is given by 



For a top-hat filter of width R we have that 
1 3ji(fci?) 



W k = 



(2tt) 3 / 2 kR 



(C5) 



where j\ is a spherical Bessel function of order 1 and is 
explicitly given by 



sin x — x cos x 



It follows that 



4 = / fr s (k) 



3h(kR) 



1 2 



kR 



(C6) 



(C7) 



2. Determination of A 
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FIG. 7: Dependence of A on the cosmological parameters, 
normalized to its value for Aao = 0.73. 



The power spectrum determined by inflation is CMB 



normalized, which fixes the value of A (see Eq. (4.6)) 



The power spectrum of the gravitational potential is 
given by Eq. (4.4) which implies that the matter power 



spectrum at redshift z = is given, after use of the Pois- 
son equation (3.9) by 



V 5 {k) = 



k 4 



y s 'm0 u /7 



(C8) 



S n( x > *) = / d3 y 6 (Vi *) Wr(\x -y\). (C2) It follows that the value of a 2 determine d from this spec- 



Decomposing the window function in Fourier modes W k , 
as in Eq. (CI), it is easily concluded that the Fourier 



trum is, using the definitions of Section IV A given by 

2 



components of the smoothed density field are 

S R (k,t) = (27r) 3 / 2 S k W k . (C3) 



~2 _ 2 
8 = a R 8 = 



dk k 4 



2^ 
5 goo 



A 



It is then clear that the variance a 2 R = S^Xft) can be 
expressed as 

4 = / T& Ps{k)w2k = I T Vs{k)wl (C4) 



k H 4 



k 

kcMB 



T 2 (k) 



3ji(fci? 8 ) 

fci? 8 



(C9) 



with i?g = 8h 1 Mpc. This value of erf depends on 
the cosmological parameters, in particular via the ratio 
go /goo and fc eq , and on n s . 
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In order to determine A, we impose that the two power 
spectra give the same erf. With the notations of Sec- 
tion [IVB] we have that 

Akt ' <l " 



2w 2 (1 



u 2)3/2 



3ji(au) 



au 



(CIO) 



with a = R 8 k c = 8ft- x Mpc x /i/20 Mpc -1 = 2/5. This 
latter integral evaluates to ~ 0.158 so that Ak\ ~ 
6.32 do. We thus use 



Such a choice allows us to tune the APM power spec- 
trum to match the inflationary power spectrum with the 
drawback that it is not CMB normalized (see Fig. [I]) . As 
explained in the text, this is not too much of a prob- 
lem since we only use it to derive analytical orders of 
magnitude and since the computations can be performed 
numerically with the inflationary power spectrum. The 
scaling of A with the cosmological parameters is depicted 
in Fig. [7} 



A = 6.32fe-V|(A,Q m o,0 AO ,n s ) 

~ 1.01 X 10 6 /l" 4 cr|(A,an0,^A0,^ S )MpC 4 (Cll) 
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